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Abstract: We analyze the possibility of studying the heavy flavour hadropro- 
duction properties at the HERA-B experiment. In addition to the high statistics 
single inclusive B spectra measurements, the measurement of the BB meson cor- 
relations is considered. The techniques of momentum estimators, widely used in 
the charm sector, are demonstrated to be useful for the BB correlation studies at 
HERA-B. The kinematic limits for the precision of the momentum estimator within 
which the pair spectra can be measured are determined. The errors are weakly 
dependent on the topology of the multibody B meson decay. 

1 Introduction 

Although the flagship goal of the HERA-B experiment is to study the CP violation in the B 
system, there should be quite an interesting data sample available for studies of production 
properties of heavy hadrons. Here we report on a study of possibilities for measurements of 
various heavy hadron production properties with an emphasis on the measurements of the 
distribution of pairs. We concentrate on beauty hadrons, as the momentum estimators were 
already demonstrated to be a useful tool in the charm sector. First, we briefly review the 
existing experimental results on the fixed target hadroproduction of heavy hadron pairs (mostly 
for the charm sector) and we briefly discuss the momentum estimator techniques used. Next, 
we find the lower boundaries for the systematic errors for the distributions of pairs when the 
momentum estimator is used in the kinematic situation of HERA-B. This is illustrated for a 
typical B decay channel. The major systematic effect on the momentum estimator coming 
from the secondary vertex resolution is also evaluated. Finally, the errors of the momentum 
estimator are studied for a wide range of B decays classified topologically. 

2 Single inclusive B production 

With a heavy hadron sample of the order of 10 5 B mesons per year the HERA-B will be 
a very competitive experiment for the study of properties of the single inclusive production. 
Although the QCD calculations of the single inclusive heavy hadron production have been 
performed some time ago [0], the next-to-leading order (NLO) QCD calculations || must be 
supplemented by various soft phenomena in order to give predictions that agree with the 



experimental data. There are two sources of theoretical uncertainties: the choice of the heavy 
quark masses, the renormalization and factorization scales etc., which are input for the NLO 
QCD calculations, as well as the parameters of the nonperturbative models which supplement 
them. Thus, the resulting total cross-section for the bottom production at HERA-B energies 
is predicted only within a range of 6-17 nb. The single inclusive Feynman-a; and the transverse 
momentum spectra are predicted with similar uncertainties 0]. 

On phenomenological grounds, the high statistics sample of B mesons may be used to answer 
questions concerning the xp and pt spectra for heavy hadrons. Namely, several experiments 
(see e.g. || and references cited therein) report non-vanishing x° F , the center of Feynman-x 
distribution 
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This offset, although predicted by the theory, is experimentally demonstrated only within one 
standard deviation. Another problem is the shape (of the tail) of the spectrum. While the 
e -b2P T f orm i s used for most of the data, some collaborations find the simple exponential form 
e -bip T better describes the data, particularly in the high pt region (see e.g. [Q). 



3 Current results BB correlations 



The perturbative and nonperturbative input parameters produce considerable uncertainties for 
the predictions of the single inclusive heavy hadron production. The distributions of heavy 
hadron pairs seems to be even more sensitive to the choice of those parameters. In addition, 
some nontrivial effects for distributions of pairs, absent in the leading order QCD (where the 
heavy hadron pair is produced back-to-back), arise entirely as the NLO or the nonperturbative 
effects. Thus, the distributions of heavy hadron pairs, although more difficult from experi- 
mental point of view, are an excellent place to study the effects of higher order QCD and 
nonperturbative contributions. 

The kinematic variables commonly used to study the hadroproduction of heavy hadron pairs 
may be divided into two groups 



• the angular variables, e.g. 

- azimuthal angle difference A</> = \<fii — 02 1, 

- pseudorapiditj difference At] = 1 771 — 772 1 ; 

• the momentum variables, e.g. 

- the Feynman-x of a pair x V p ir , 

- the effective mass of the pair M e ff, 

- the transverse momentum of the pair p^ imT = (pit + P2t) 2 , 

- the rapidity difference y<uff = \yi — 2/2 1 • 



A majority of data comes from the fixed target hadroproduction of charmed pairs. The 
only result for bottom is based on 9 pairs (not fully reconstructed). The results are collected 
in Table 1. In all these studies momentum estimators have been used. 
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Table 1: Experimental results for the heavy hadron momentum correlations. The first four are 
the results for the charm pairs, the last one is for the beauty pairs. 



4 Momentum estimators 



To study the distributions of angular variables of the heavy hadrons one needs to know the 
flight directions of both hadrons. This is achieved from the precisely determined primary and 
secondary vertices. It is enough that one (or even none) of the hadrons is fully reconstructed. 

In order to study the momentum correlations, one needs to know the momenta of both 
heavy hadronsQ. There are various techniques used for the estimation of the momentum of 
decaying heavy hadrons. They are based on the precise measurement of the position of the 
primary and secondary vertex and the measurement of the charged decay products' momenta. 
The kinematic techniques described below are often supplemented by Monte Carlo information. 
There are two popular approaches: 



• estimator EQ, the problem may be solved exactly with three inputs: the heavy hadron 
flight vector, the visible decay products' momenta and the effective mass of the invisible 
decay products. 

• estimator ET, relies on the assumption that the invisible decay products' momentum, 
in the rest frame of the heavy hadron, is perpendicular to its flight vector. Then the 
corresponding boost from the rest frame may be found to match the laboratory visible 
momentum. 



The estimator EQ gives the exact answer, which is ambiguous. Then either the most probable 
solution (as suggested by MC or the ET estimator) or just the average of the two solutions 
may be used. For this estimator one needs to know the effective mass of the invisible decay 
products. This is simple when the invisible decay product is e.g. one 7r°, neutrino or a narrow 
state decaying into neutrals. But in a case of e.g. nonresonant ttV in the final state, the 
necessary estimates of the invisible effective mass and the experimental uncertainties of the 

1 There is an alternative approach to study the correlations of heavy hadron pairs in a case of experiment 



like HERA-B 12 . In this report we concentrate on the momentum estimator approach. 



measured quantities needed for the exact solution, limit the applicability of this estimator to a 
small class of events. 



The ET estimator is simple to apply and works well for multibody final states Q. It has 
been used for studies of the single inclusive distributions by the E653 collaboration [13|, and 
the charm hadron momentum correlations in the E653 ||, WA75 0, WA92 || and NA32 || 
data. 



5 HERA-B and BB correlations 

According to expectations fl|], the HERA-B experiment should collect a sample of about 10 5 
BB events per year. The advantage of the design of the experiment is that one of the B 
hadrons in each event will be fully reconstructed. The other B hadron will decay into any 
decay mode. The momentum estimator ET has typically a sizable error, however, the es- 
timated momentum of the other B will be added to the well measured, fully reconstructed 
B. Thus, the relative error of the momentum of the BB pair will be roughly halved. The 



results of the simulation [f[4|, with HERA-B parameters and xf and pr acceptances 0, for 
the BB pairs with one of the B mesons fully reconstructed and the other decaying into 
B~ — > .D* + (2010)7r - 7r~7r — > (i-C~7r + 7r )7r + 7r~7r~7r chain (D° decay products are in the brackets 
and both 7r°s are unmeasured) are shown in Fig.l. The resolution of the momentum estimator 
for the laboratory momentum of the single B meson is 18 %, while for the laboratory momen- 
tum of the pair pf^ r is 9.1%. This results in the errors for the pair variables^ Ax^ air =0.027, 
A^ mT= i. 3 GeV 2 , AM eff = 0.34 GeV and Ay diff =0.1Q. 



Figure 1: The relative difference of the estimated and the true laboratory momentum for the 
single B meson (left) and for the BB pair (right). One of the B mesons in the pair decays into 
B~ — ► .D* + (2010)7r~7r - 7r — ► (K~7r + 7r°)7r + 7r~7r~7r with both n°s being unmeasured. The other 
B meson is fully reconstructed. 



The quoted errors of the pair variables are the lower limits, determined by the produc- 
tion and the decay kinematics, when the momentum estimator ET is used. One of the most 
important sources of additional systematics is the secondary vertex resolution^. To use the 
momentum estimator, only the flight vector direction, and not e.g. the decay length, must be 
known. The longitudinal resolution has a little influence on the flight vector direction. This is 

2 In practical applications, in order to reduce the error of the estimator, the measured distributions are cut at 
the very tail (with a typical loss of a few percent of events) . A strong improvement for the momentum estimator 
comes from the application of the cut on the transverse momentum of the visible decay products with respect 
to the flight vector. This cut, however, has much stronger influence on the statistics. For the following studies 
of the errors of the momentum estimator we do not optimize the selection cuts, thus keeping more than 90% of 
events. 

3 The B mesons in the BB pair are uncorrelated, thus giving the averages < p s ^ m T >=10.3 GeV 2 and 
< M eff >=11.8 GeV. 

4 The primary vertex resolution at HERA-B is of the order of 10/im jij and is neglected here. 
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Table 2: Errors of the momentum estimator ET for various classes of the B decaying into the 
+ (wr) channels. In brackets are the decay products of the D^*\ The Tt s are not being 
measured. The other B meson in the event is fully reconstructed. 



further reduced if the decay length cut is applied. Indeed, the simulation with the secondary 
vertex resolution a z = 500/im and a x = a y = 25/im JIJ results in small changes of the momen- 
tum estimator errors for the pair Feynman-x, Ax P p tr =0.028, and for the rapidity difference, 
A?/di//=0.17, once the cutf] of 3 mm on the decay length is applied. The effect is stronger for 
the transverse momentum Ap s 2 mT =1.8 GeV 2 and the effective mass AM e ff = 0.38 GeV. 



6 Topological approach 

As shown in the previous section, the ET momentum estimator is adequate within a typical 
error of ~10% for the laboratory momentum of the BB pair. The studied decay channel has 
the branching ratio of (1.5±0.7)% |H| only. On the other hand, the inclusive branching ratios 
for B mesons are B -> D+X = (26 ± 4)% and B -> D°X = (54 ± 6)% pf. As the B mesons 



decay in a large number of channels with small branching ratios, in any practical approach, 
one needs to identify the decay vertex of the other B hadron and all tracks which belong to 
this vertex. Exact knowledge of the decay channel and the decay chain is not necessary for the 
momentum estimator ET to be applied. The dependence on the decay chain for the decays 
with the same number of tracks emerging from the secondary vertex is week (see Table 2). In 
addition, the vertex missing mass may be used to estimate the number of missing neutrals. 
The events may be classified topologically and their systematics studied for the whole class of 
events. 



7 Conclusions 

We show that the HERA-B experiment will provide an excellent opportunity to study the heavy 
quark hadroproduction and the perturbative and nonperturbative contributions to this process. 
It should resolve long standing problems of the single inclusive heavy hadron production like 
the shape of the Feynman-x and px distributions. The heavy hadron pair distributions are 
also worth studying at HERA-B. Firstly, there should be a large BB sample for the study 

5 The mean flight path of the B meson at HERA-B energies is of the order of 9 mm. 



of angular correlations (comparing with current results). In addition, the study of the BB 
momentum correlations, although difficult, is possible with the use of momentum estimators. 
Approximately a few percent of BB events, with the first B meson fully reconstructed and the 
second B meson partially reconstructed, could be used, compared with <C 1CT 3 , when both B 
mesons are required to be fully reconstructed. The kinematic limits on the resolutions of the 
momentum estimator are of the order of 8-12% for the laboratory momentum of the BB pair, 
resulting in the resolutions 

• 0.02-0.04 for the x P p ir distribution, 

• 1.2-1.4 GeV 2 for the p s l mT distribution (with the average < p s l mT >=10.3 GeV 2 ), 

• 0.3-0.5 GeV for the M eff distribution (with the average < M eff >=11.8 GeV), 

• 0.15-0.20 for the r/diff distribution. 

The source of the main systematic error, namely the effects of the secondary vertex resolution, 
are under control once the cut on the secondary vertex separation of 2-3 mm is applied. 
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